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Destabilization of rotating flows with positive shear by azimuthal magnetic fields 
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According to Rayleigh’s criterion, rotating flows are linearly stable when their specific angular 
momentum increases radially outward. The celebrated magnetorotational instability opens a way 
to destabilize those flows, as long as the angular velocity is decreasing outward. Using a local 
approximation we demonstrate that even flows with very steep positive shear can be destabilized 
by azimuthal magnetic fields which are current-free within the fluid. We illustrate the transition of 
this instability to a rotationally enhanced kink-type instability in case of a homogeneous current in 
the fluid, and discuss the prospects for observing it in a magnetized Taylor-Couette flow. 
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From the purely hydrodynamic point of view, rotat¬ 
ing flows are stable as long as their angular momentum 
is increasing radially outward [l]. Since this criterion 
applies to the Keplerian rotation profiles which are typ¬ 
ical for low-mass accretion disks, the growth mechanism 
of central objects, such as protostars and black holes, 
had been a conundrum for many decades. Nowadays, 
the magnetorotational instability (MRI) [ 2 ] is considered 
the main candidate to explain turbulence and enhanced 
angular momentum in accretion disks. The standard 
version of MRI (SMRI), with a vertical magnetic field 
B z applied to the rotating flow, requires both the rota¬ 
tion period and the Alfven crossing time to be shorter 
than the timescale for magnetic diffusion [3j. This im¬ 
plies, for a disk of height H , that both the magnetic 
Reynolds number Rrn = poaH 2 Lt and the Lundquist 
number S = hquHva must be larger than one (fl is the 
angular velocity, p ,0 is the magnetic permeability con¬ 
stant, a the conductivity, va '■= B z /y/poP is the Alfven 
velocity, with p denoting the density). While these condi¬ 
tions are safely fulfilled in well-conducting parts of accre¬ 
tion disks, the situation is less clear in the “dead zones” of 
protoplanetary disks, in stellar interiors and liquid cores 
of planets, because of the small value of the magnetic 
Prandtl number Pm = v/p Q, i.e. the ratio of viscosity 
v to magnetic diffusivity rj := (/j, 0 ct) _1 . 

This low Pm case is also the subject of intense theo¬ 
retical and experimental research initiated by Hollerbacli 
and Rudiger :5|. Adding an azimuthal magnetic field B^ 
to B z , the authors found a new version of MRI, now 
called helical MRI (HMRI). It was proved to work also 
in the inductionless limit [6j, Pm = 0, and to be gov¬ 
erned by the Reynolds number Re = RmPm - and the 
Hartmann number Ha = 5Pm^^ 2 , quite in contrast to 
standard SMRI that is governed by Rm and S. 

A somewhat soberinglimitation of HMRI was iden¬ 
tified by Liu et al. [7| who used a local approxi¬ 
mation (also called short-wavelength, Wentzel-Kramers- 


Brillouin (WKB), or geometric optics approximation, 
see Q) to find a minimum steepness of the rotation 
profile n(r), expressed by the Rossby number Ro := 
r(2H)- 1 9H/ l 9r, of Ro L ll = 2(1—v/2) « -0.828. This 
lower Liu limit (LLL) implies that, at least for B<f,(r) oc 
1/r, HMRI does not extend to the most relevant Keple¬ 
rian case, characterized by Roxep = —3/4. Surprisingly, 
in addition to the LLL, the authors found also a second 
threshold of Ro, which we call upper Liu limit (ULL), 
at Roull = 2(l+\/2) ~ +4.828. For Ro > Roull one 
expects a magnetic destabilization of those flows with 
strongly increasing angular velocity that would even be 
stable with respect to SMRI. 

By relaxing the demand that the azimuthal field is 
current-free in the liquid, i.e. B^(r) oc 1/r, and allowing 
fields with arbitrary radial dependence, we have recently 
shown [s, [i)f that the LLL and the ULL are just the end¬ 
points of one common instability curve in a plane that is 
spanned by Ro and a corresponding steepness of the az¬ 
imuthal magnetic field, called magnetic Rossby number, 
Rb := r(2B ( / > /r)~ 1 d(B ! j > /r)/dr. In the limit of large Re 
and Ha, this curve acquires the closed and simple form 


Rb = - 


1 (Ro + 2) 2 
8 Ro + 1 


( 1 ) 


A non-axisymmetric “relative” of HMRI, the az¬ 
imuthal MRI (AMRI) 0, which appears for purely or 
dominantly B has been shown to be governed by basi¬ 
cally the same scaling behaviour, and the same Liu limits 
EL Actually, the key parameter dependencies of HMRI 
and AMRI were confirmed in various liquid metal exper¬ 
iments at the PROMISE facility 12, Kl¬ 
in the present paper, we focus exclusively on the case 
of positive Ro, i.e. on flows whose angular velocity (not 
only the angular frequency) is increasing outward. From 
the purely hydrodynamic point of view, such flows are 
linearly stable (while non-linear instabilities were actu¬ 
ally observed in experiments 14]). Flows with positive 
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Ro are indeed relevant for the equator-near strip (ap¬ 
proximately between ±30°) of the solar tacliocline [iff, 
which is, interestingly, also the region of sunspot activ¬ 
ity [ill- Up to present, the ULL at Roull = +4.828 
has only been predicted in the framework of various lo¬ 
cal approximations while attempts to confirm it 

in a 1-dinrensional modal stability code on the basis of 
Taylor-Couette (TC) flows have failed so far 0 . Hence, 
the questions arise: Is the magnetically triggered flow 
instability for Ro > Roull a real phenomenon (which 


would fundamentally modify the stability criteria for ro¬ 
tating flows in general), or just an artifact of the local 
approximation, and is there any chance to observe it in 
a TC experiment? 

In order to tackle these problems we restrict our atten¬ 
tion here to non-axisymmetric instabilities, which are the 
relevant ones for pure B and further assume Pm = 0. 
Under these assumptions, we had recently [§] derived the 
closed equation 


1 [(1 + HaV) 2 - 4Ha 2 Rb(l + HaV) - 4HaV][l + Ha 2 (n 2 - 2Rb)] 2 
4 Ha 4 RoV - [(1 + Ha 2 (n 2 - 2Rb)) 2 - 4HaV][Ro + 1] 


for the marginal curves of the instability, where the fol¬ 
lowing definitions for Re, Ha and the modified azimuthal 
wavenumber n are used: 


Re = 


Ha = 


a f l(r) 


|k | 2 

a 




|k | 2 r(poppv ) 1 / 2 

n = m/a , 


( 3 ) 

( 4 ) 

( 5 ) 


with a = fc 2 /|k| and |k | 2 = fc 2 + fc 2 defined as functions 
of the axial and radial wavenumbers k r and k z . 

Because of its comparably simple form, and the ab¬ 
sence of the ratio /? of azimuthal to axial magnetic field 
(which would play a decisive role for HMRI), Equa¬ 
tion (2) allows to easily visualize the transition from a 
shear-driven instability of the AMRI-type to the current- 
driven, kink-type Tayler instability (TI) jl8} , when going 
over from Rb = — 1 to Rb = 0. 

Let us start with the current-free case, Rb = — 1. Fig¬ 
ure la shows, for varying values of Ro and the particular 
case n = 1.4, the marginal curves in the Ha-Re plane. We 
see that the critical Re increases steeply for Ro below 6 
which reflects the fact that we approach Roull = 4.828 
from above. We ask now for the dominant wavenum¬ 
bers, as illustrated in Figure lb for the particular value 
Ro = 5.5. Evidently, the minimal values of Re and Ha 
(the “knee” of the curve) appear for n ~ 1.4 which 
represents a rather “benign” combination of wavenum¬ 
bers with k r ~ k Zl so that neither the axial nor the ra¬ 
dial wavelength of the perturbations diverges. From this 
point of view, there seems to be no contradiction with 
the underlying short-wavelength approximation. 

While for Rb = — 1 the only energy source of the in¬ 
stability is the shear of the rotating flow, we move now in 
direction of Rb = 0 which corresponds to a constant cur¬ 


rent density in the fluid, for which the kink-type TI 18] 


is expected to occur. For the particular choice n = 1.2, 


this transition is illustrated in Figure 2 where we have 
intensionally chosen, for all Rb, the same scales for Re 
and Ha. For Rb = —0.6 we observe the appearance of 
a crossing with the abscissa, i.e. a point where the in¬ 
stability draws all its energy from the electrical current 
instead of the shear. Actually, the lowest value where 
this can occur is Rb = n 2 /4 — 1 = —0.64 

For Rb = 0 the instability is characterized in more 
detail in Figure 3. Very similar to the results of [Hj], we 
observe in Figure 3a that for Ro > 0 the curves move 
to the left with increasing Re (i.e., the flow supports the 
kink-type instability) and converge to well-defined values 
of Ha when Re goes to infinity. The dependence on the 
wavenumber ratio a is quite interesting. Figure 3b shows 
that the mode with n = 1 (i.e. with k r = 0), which is still 
dominant at Re = 0, is replaced by modes with higher 
values of n for increasing Re. The limits of the critical Ha 
for Re = 0 and Re —)• oo can be determined by setting 
to zero, in Equation 2, the nominator or denominator, 
respectively, which leads (for Rb = 0) to 


Ha Re=0 = \/\Jn{2 - n) 

HaRe-foo — 


(Ro + 1) + \J (Ro + l)(Ro + 2 )/n 


Ro + (Ro + 1)(4 — n 2 ) 


( 6 ) 

■( 7 ) 


In the limit Ro —> oo the limit values of Ha converge 
slowly to zero according to Ha( RejRo )_,. 00 ~ n _1 / 2 Ro -1 ^ 4 . 

In the following, we compare our WKB results with re¬ 
cent findings [h| obtained for a TC flow with inner and 
outer radii r, and r Q rotating with the angular velocities 
fl; and S2 0 , respectively. The corresponding ratios are de¬ 
fined as fj = ri/r 01 and (i = f2 0 /fV For this TC configu¬ 
ration, the following modified definitions of the Reynolds 
and Hartmann number were used: Re = H 0 ri(r 0 — Ti)/v , 
Ha = H^(ri)(ri(r 0 - ri)) 1 / 2 /(popup) 1 / 2 . The non-trivial 
point is now how to translate the p of a TC flow, char¬ 
acterized by H(r) = a + b/r 2 , to the Ro of a flow with 
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FIG. 1: Marginal curves for Rb = —1. (a) Dependence on Ro 
for n = 1.4. (b) Dependence on n for Ro = 5.5. The insert 
shows the dependence of the minimum value (with respect to 
Ha) of the critical Re on n. The arrow points to the optimum 
n « 1.35 that leads to the lowest critical Re. 


fi(r) ~ r 2Ro . An often used correspondence, based on 
equalizing the corresponding angular velocities at rj and 
r 0 fill • leads to 


Ro* ~ —l/21ogrj/t (8) 


while an alternative, more shear-oriented version leads to 


* 1 (1 + 77) (A - 1 ) 

“ 2 ( 1 -»))(£ + !) ' 


(9) 


Actually, for comparably small (positive or negative) val¬ 
ues of Ro, the differences are not very significant, but 
they increase for steeper profiles. This is a key point for 
the adequateness of TC flows to “emulate” steep power 
function flows. In 0 , the destabilizing effect of positive 
shear had been studied for TC flows (with Rb = 0 only), 
both for a wide gap with r) = 0.5 as well as a narrow 
gap with 77 = 0.95. In either case, for large values of 
/t, the critical Ha converged to some non-zero constant, 
which is not compatible with the translation to Ro* since 
the latter should lead to a zero critical Ha (according to 


0 1 


Rb=-1.0 


0 12 3 4 

Rb=-0.: 


0 12 3 4 


0 12 3 4 

Rb=-0.3 




ya 

1 ■ 
0 

Ci 





* 



\ 



_ 





FIG. 2: Marginal curves for n = 1.2 and various values of Rb, 
in dependence on Ro. From top to bottom, the instability 
changes its character from a (magnetically triggered) shear- 
driven instability to a (rotationally influenced) current-driven 
TI. For n = 1.2, TI appears first for Rb = n 2 /4 — 1 = —0.64. 
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FIG. 3: Marginal curves for Rb = 0. (a) Dependence on Ro 
for n = 1.0. (b) Dependence on n for Ro = 5. 


H a (Re,R 0 )->oo — n 1//2 Ro 4 ' 4 , see above). It turns out 
that the translation to Ro** is physically more adequate. 

With the reasonable choice k z = k r = n /(r* — r 0 ) we 
obtain the translations Re = tt 2 2 5 / 2 fifj/((1+— Re 
and Ha = 7r 2 (l + fj) 2 / {{2fjY^(l — i)) 3 / 2 ) Ha. For 
fj = 0.95 this amounts to Re = 1061/(1 + l//t)Re and 
Ha = 2435Ha. Figure 4 shows the corresponding WKB 
results, both for assuming a translation to Ro* (dashed 
lines) and to Ro** (full lines). For Re = 0 our result 
Ha = 2670 agrees reasonably well with the exact value 
Ha = 3060 of the modal stability analysis []jj. What is 
more, the typical bend of the marginal curve to the left 
for increasing Re, and the limit values of Ha for large Re, 
are also confirmed. Yet, subtle differences show up for 
the two ways of translation: the use of Ro** confirms the 
existence of a finite limit value for the critical Ha, as typ¬ 
ical for TC flows, while the use of Ro* would ultimately 
lead to a zero limit value. 

This encouraging consistency of the local approxima¬ 
tion and the modal stability analysis, evidenced for Rb = 
0, brings us back to the point whether, for Rb = — 1, the 
ULL can be confirmed in a TC experiment. Assuming 



Ha 

FIG. 4: Mar gina l curve for Rb = 0 and n = 1.41, as scaled 
according to [l9j ]. The full lines correspond to the translation 
of jj, to Ro**, the dashed lines to Ro*. 


Ro** as more physical than Ro*, in the limit /t —> oo we 
obtain Ro-1^ = 1/2(1-H))/(1 — fj). This means, in turn, 
that to emulate some Ro in a TC-flow, fj has to fulfill the 
relation fj = (2Ro — l)/(2Ro + 1). With view on the 
ULL, this implies that for Ro = 6, say, a minimum value 
of fj = 11/13 = 0.846 is needed. For TC-flows with wider 
gaps, such as fj = 1/2 the necessary shear could simply 
not be realized. 

What are, then, the prospects for a corresponding 
experiment? Evidently, we need a rather narrow gap 
flow. Let us stick, for a first estimate, to the safe value 
fj = 0.95, and take the typical values Ro = 6, Ha = 2 and 
Re = 12 as read off from Figure la. This translates to 
/t = 1.89, Re = 8324 and Ha = 4870. For a prospective 
TC experiment with Na at 150°C, with p = 910 kg/m 3 , 
v = 5.94 x 10 -7 m 2 /s, a = 9 x 10 6 S/m, and an outer 
diameter of r 0 = 0.25 m, this would amount to a rather 
moderate rotation frequency of fl D / ( 27 t) = 0.26 Hz, yet a 
huge magnetic field = 0.69 T that requires a cen¬ 

tral current of I = 8.6 x 10 5 A. Exhausting the shear re¬ 
sources, by choosing jl —> oo and fj = 0.85 « 11/13, those 
values would drop to Re = 3796, Ha = 892 or, physically, 
to n o /( 2tt) = 0.044 Hz, B 4> {n) = 77 mT, I = 8.2 x 10 4 
A. Any real TC experiment, however, would need more 
detailed simulations with a ID marginal stability code to 
confirm and optimize the parameters. 
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